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基于 火焰 面 /进度 变量 方法 模拟 超 临界 
压力 下 的 液 氧 甲 迷 火焰 
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摘 要 :为 了 深入 认识 超 临 界 压 力 下 低温 推进 剂 的 复杂 燃烧 过 程 ,基于 OpenFOAM 开源 CFD 仿真 
ds A es be 一 套 超 临界 We 
方法 。 针 对 Mascotte 试验 台 的 液 氧 甲烷 燃烧 实验 进行 了 数值 模拟 ,结果 较 准 确 地 捕 提 了 火焰 形态 
同时 研究 了 燃烧 室 压 力 对 火焰 的 影响 。 研 究 表明 :由 于 超 临 界 压力 下 液 所 在 升温 时 存在 拟 沸 腾 过 
) 程 ,射流 剪 切 层 中 液 氧 迅速 膨胀 并 产生 了 径 向 向 外 的 流动 ,造成 了 火焰 的 突然 扩张 并 在 火焰 下 游 形 
成 回流 区 ;燃烧 室 压 力 升 高 时 甲烷 密度 增加 而 入 口 速 度 减 小 ,同时 液 氧 升 温 膨 胀 的 幅度 减 小 ,造成 
了 火焰 长 度 的 增加 以 及 火焰 下 游 回流 区 的 消失 。 
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Simulation of LO, /CH, flame under supercritical pressure 
based on flamelet/ progress variable method 


LI Yuhang , WANG Yi' ,WU Baoyuan’ , YANG Jianwen' 


(1. Key Laboratory for Liquid Rocket Engine Technology, Xi’an Aerospace Propulsion Institute ,710100 Xi’an, China; 
2. Academy of Aerospace Propulsion Technology ,710100 Xi’an , China) 


Abstract : In order to get an in-depth insight into the complicated combustion process of cryogenic propel- 
lant under supercritical pressure ,a homogeneous simulation method of supercritical combustion is devel- 
oped using open source CFD toolbox OpenFOAM, based on real-fluid thermo-physical model and flamelet/ 
progress variable combustion model. The LO,/CH, flame of Mascotte combustor is simulated, and the effect 
of chamber pressure on flame structure is studied. Results show that due to the existence of pseudo-boiling 
process during LO, heating up under supercritical pressure , the LO, in shear layer expands rapidly along 
radial direction, leading to the reverse flow region downstream of the flame. During the rise of chamber 
pressure ,the density of methane increases with its inlet velocity decrease ,and the expanding extant of LO, 
reduces ,thus resulting in the increase of flame length and the vanishing of reverse flow region downstream 
of the flame. 
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第 5 期 EGR, E: 基于 火焰 面 /进度 变量 
大 推力 液体 火箭 发 动机 的 主 燃烧 室 压 力 可 达 
20 MPa 以 上 , 远 高 于 液 氧 .甲烷 .煤油 等 典型 推进 剂 
的 临界 压力 。 在 超 临 界 压力 条 件 下 ,流体 的 表面 张 
力 及 气 液 相 界面 消失 ,低压 下 经 典 的 液体 雾 化 破碎 
过 程 转 变 为 一 种 由 扩散 过 程控 制 的 稠密 流体 的 挨 混 
行为 。 这 样 的 超 临 界 流动 现象 已 被 一 系列 实验 研究 
所 证 实 。 德 国 宇航 中 心 (DLR ) 与 美国 空军 研究 实验 
室 (AFRL) 先 后 对 液 氮 单 射流 、 液 所 /氯气 同 轴 射 流 
进行 了 系统 的 实验 研究 1 ,观察 到 了 高 压 下 液 滴 破 
碎 现 象 的 消失 ,并 对 射流 长 度 .扩张 角 以 及 密度 分 布 
进行 定量 测量 。 法 国 EM2C 实验 室 建 立 了 Mascotte 
试验 台中 ,针对 LOVCH4 与 L0H, 超 临界 燃烧 进行 
光学 测量 ,给 出 了 火焰 的 -OH 图 像 。 
二 航 空 航天 工程 中 开展 试验 的 成 本 高 昂 , 通 过 计 
算 沛 体力 学 方法 进行 研究 和 辅助 设计 可 以 降低 时 间 
PARRAK, 且 能 够 提供 更 详尽 而 全 面 的 数据 。 
已 前 对 于 发 动机 相关 条 件 下 的 流动 与 燃烧 数值 模拟 
EE 初步 形成 体系 , 即 采用 真实 流体 模型 的 均 相 方 
PS) ,但 对 于 燃烧 过 程 的 处 理 仍 面临 许多 挑战 。 一 
充 面 ,真实 流体 状态 方程 及 化 学 反应 计算 的 非 线性 
特 任 带 来 了 程序 收敛 性 方面 的 困难 ; 另 一 方面 ,多 步 
化 村 反 应 机 理 涉及 的 刚性 常 微分 方程 组 求解 使 得 化 
学 及 应 源 项 计算 过 程 非常 缓慢 。 目 前 大 多 数 研究 都 
避 碗 采用 有 限 速率 反应 模型 结合 详细 化 学 反应 机 理 
进 得 仿真 。 文 献 [11] 采 用 6 组 分 7 步 化 学 反应 机 理 
研 容 了 LO/H, 火 焰 。 文 献 [12] 采 用 单 步 总 包 化 学 
反 庶 机 理 研究 了 LO./CH, 火 焰 , 并 使 用 了 修正 方法 
弥 神 总 包 反 应 在 火焰 温度 及 层 流 火焰 传播 速度 预测 
上 的 不 准确 。 文 献 [13-14] 采 用 无 限 快速 化 学 反应 
模型 分 别 研究 了 单 喷嘴 和 多 喷嘴 的 LO / CH, RABE 
程 。 文 献 [15] 采 用 扩展 平衡 态 方法 研究 了 42 喷嘴 
火箭 发 动机 的 燃烧 , 较 好 地 预测 了 推力 室 中 的 燃烧 
不 稳定 。 然 而 ,无 论 是 简化 反应 机 理 或 无 限 快速 反 
应 模型 ,都 对 化 学 反应 计算 引入 了 过 多 的 近似 。 
火焰 面 模型 作为 一 种 新 型 的 燃烧 建 模 思路 ,能 
够 在 考虑 详细 化 学 反应 机 理 的 同时 ,通过 预先 建 表 
过 程 大 幅 减 小 计算 量 。Peters 最 早 提 出 稳 态 层 流 火 
焰 面 (steady laminar flamelet, SLF ) HEA! , 并 讨论 
了 该 模型 所 适用 的 燃烧 模式 。Pierce 于 2004 年 做 出 
里 程 碑 式 的 改进 ,使 用 反应 进度 变量 描述 化 学 反应 
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(flamelet/progress variable, FPV) EW] 。 此 方法 具 
有 捕 换 火焰 抬 举 局 部 熄火 等 非 稳 态 人 燃烧 现象 的 能 
力 , 目 前 已 在 航空 发 动机 、 超 燃 冲压 发 动机 仿真 方面 
得 到 广泛 应 用 。 近 年 来 ,许多 研究 者 开始 发 展 火焰 
面 模型 在 火 第 发 动机 相关 超 临 界 燃烧 仿真 方面 的 能 
力 。 文献 [18] 采用 稳 态 火焰 面 模 型 研究 了 LO,/ 
CHs 超 临界 同 轴 射 流 火焰 。 文 献 119-20 ] 阐述 了 火 
焰 面 /进度 变量 模型 与 大 涡 模 拟 相 结合 的 超 临 界 燃 
烧 仿 真 方法 体系 。 文 献 [21-23 ] 采 用 火焰 面 模型 对 
液 氧 煤油 发 动机 所 使 用 的 气 液 同 轴 离 心 式 喷 嘴 E 
轴 双 离心 式 喷嘴 的 燃烧 过 程 进行 了 一 系列 仿真 工 
VE ,对比 了 关键 几何 参数 对 于 燃烧 特性 的 影响 规律 ， 
得 出 了 许多 具有 工程 价值 的 结论 。 

本 研究 将 真实 流体 热 物 性 植 和 人 火焰 面 / 进 度 变 
量 方法 ,基于 OpenFOAM 开源 仿真 平台 开发 了 低温 
推进 剂 超 临 界 燃烧 的 均 相 求解 器 。 针 对 Mascotte 燃 
烧 试 验 台 的 液 氧 甲 烧 超 临界 火焰 进行 了 数值 模拟 ， 
验证 了 方法 的 准确 性 ,并 研究 了 燃烧 室 压 力 对 于 火 


焰 的 影响 。 


1 数学 物理 模型 


1.1 流 场 控制 方程 


本 研究 基于 RANS 方法 和 FPV 模型 进行 燃烧 模 
拟 , 控 制 方程 涉及 到 Favre 平均 的 质量 、 动 量 、 混 合 
分 数 .混合 分 数 方差 及 反应 进度 变量 的 输 运 方程 。 


B+v. (pic) =0 a) 
t + (puu) = —Vp+¥ 5 (2) 
My. pif) =v [ (a0 +5) (3) 


apf” 3) By 
rad - (pif) = 


y- [ (PD a] +2 aw Vi-py (4) 


we ey -(puC) =V - [到 + 全 j75] + pi 
Oc 


(5) 


中 局 部 熄火 的 中 间 状 态 ,形成 了 火焰 面 进度 变量 


式 中 :D 为 分 子 扩散 系 数 ;or 与 ee 都 取 0.7;X = 
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C, TG =2. 0; im RIPE yu, 由 潮流 模型 决定 ; 进 
度 变 量 源 项 wv 从 庙 流 火焰 面 数据 库 中 查 表 得 到 。 
1.2 真实 流体 热 物 性 模型 


立方 型 状态 方程 能 够 准确 描述 真实 气体 效应 
保持 较 小 的 计算 量 。 本 研究 使 用 SRK 0 


_ RT aa(l7) 
PEV b VV +b) a 
对 于 多 组 分 混合 物 ,需要 使 用 混合 规则 计算 式 
中 的 参数 , 即 
aa = > DE (7) 
b= Yad, (8) 


0 aja; = ./a;,a,a;,0;(1 —k,) (9) 
其 中 :a,b, 分 别 为 组 分 i 的 状态 方程 参数 ; 启 为 二 元 
ADEE BE, 分 别 使 用 下 列表 达 式 计算 。 


Ei RT. 

m a, =0. 427 4g Rta) (10) 
BE Pa 

A RT, 

GN b; =0. 086 64 fa). (11) 
Nla T) = [1 + (0. 485 08 +1. 551 7w; - 

> 0. 156 13%) (1 - T%5/T%5)]? (12) 
p< AVN 了 

© k; m1- [yp (13) 


JOR, 民 为 通用 气体 常数 ;7 Pa Va w; 分 别 为 组 分 

i 的 临界 温度 ,临界 压力 ,临界 体积 以 及 偏心 因子 。 
热力 学 性 质 通过 理想 气体 参照 值 加 上 偏离 函数 

得 到 。 其 中 内 能 与 定 压 比 热 的 偏离 函数 形式 如 下 。 


(Tp) = e(n) + [LE + (e) Je (14) 


CO) = CoD) ~ [[(Sh) Je (15) 


FAP Sd BOHR GE SRK 状态 方程 进行 推导 得 到 。 
和 Co 表示 理想 气体 参照 值 ,使 用 NASA STRN 
计算 。 
APE RA Chung 方法 计算 ,其 表达 式 为 
; 1]  36.344(M + T)” 
(16) 
各 项 参数 的 计算 方法 详细 可 参照 文献 [24] 。 
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1.3 火焰 面 方程 


对 于 化 学 反应 时 间 尺 度 小 于 淇 流 时 间 尺 度 的 非 
预 混 燃 烧 过 程 , 潮 流 火焰 采用 一 维 对 撞 层 流 火焰 结 
构 的 系 综 平均 来 表征 。 采 用 Crocco 坐标 变换 可 将 
对 撞 火 焰 从 物理 空间 转换 到 混合 分 数 空间 ,得 到 包 
含 组 分 方程 与 能 量 方程 的 火焰 面 方程 组 。 

组 分 方程 形式 如 下 。 
oY; my FY; Y, 
ot ~ 2Le, af? 
RIP y =2D( af/ax)* Arp at REALS , REDRE im 
流 涡 对 火焰 的 拉 伸 作用 ,对 火焰 结构 产生 影响 ,x 为 
物理 空间 上 的 坐标 ,D 为 混合 分 数 的 扩散 系数 ;Le， 
表示 组 分 相应 的 刘易斯 数 ( 本 研究 假设 对 所 有 组 
分 均 为 1) ;w; 为 组 分 的 生成 速率 。 

能 量 方程 形式 如 下 。 

pe; 三 E af? + $ ho, + a 
式 中 :e, WEEE h, 为 组 分 对 应 的 化 学 生成 
KA o 本 研究 求解 稳 态 火焰 面 方程 , 因 因而 忽略 温度 与 
压力 的 瞬 态 项 。 


1.4 淇 流 火焰 数据 库 
注 流 化 学 交互 作用 使 用 概率 密度 函数 处 理 。 对 


pe + pw; (17) 


(18) 


前 一 节 的 层 流 火 焰 面 数据 进行 概率 密度 积分 得 到 滑 
流 火焰 面 数据 库 。 
5 = [eF O PG, afac (19) 


其 中 p KRMAR .组 分 及 热 物 性 值 。 假 
设 混 合 分 数 /与 进度 变量 C 的 概率 密度 相互 独立 ， 
则 式 中 的 概率 密度 函数 可 写 为 


P(f,C)= P(f) P(C) (20) 


假定 混合 分 数 / 服从 B 分布, 进度 变量 C 服从 5 
分 布 。 可 以 得 到 应 流 火 焰 面 数据 库 。 

ø =g(f f”, Č) (21) 

流 场 中 的 温度 、 组 分 等 数据 可 通过 混合 分 数 均 


值 和 方差 F” 以 及 进度 变量 均值 Č 等 3 个 参数 进 
行 插值 查 取 。 

图 1 给 出 了 本 人 研究 所 建立 求解 器 的 计算 过 程 ， 
其 中 真实 气体 火焰 面 方程 的 求解 与 淇 流 火 焰 面 数据 
库 的 建立 先 于 CFD 计算 而 进行 。 在 CFD 计算 时 对 
温度 .组 分 .密度 .可 压缩 性 等 数据 实时 进行 查 表 。 
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Flamelet Generation 
Solve real-gas flamelet equations 
Tool:FlameMaster 
PAS" O 
CFD Code 
Pressure-based PIMPLE algorithm 
Tool:OpenFOAM 


KIL 火焰 面 模型 的 数值 计算 流程 


Fig.1 Calculation flowchart of the present flamelet model 
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15 计算 域 模型 


oa 
会 依照 Singla 等 中 在 Mascotte 试验 台 上 所 进行 的 
L9g5CH, 实 验 设置 计算 条 件 。 该 试验 台 的 喷嘴 尺寸 

烧 室 结构 如 图 2 所 示 , 其 中 di =3.6 mm,d, = 
ie =5.6 mm,d, =10 mm, 按照 等 流通 面积 
将 原 实验 宽度 50 mm 的 方形 燃烧 室 折算 为 半径 
2821 mm 的 圆 形 燃烧 室 , 沿 切 向 截取 30° RUE WAS 
ge ee 
3 ( 琶 示 ,网 格 在 喷嘴 附近 进行 加 密 。 喷 注 器 与 燃烧 
罕 评 采用 绝热 无 滑 移 壁 面 ,出 口 采用 压力 出 口 边界 。 
表 省 列 出 了 燃料 与 氧化 剂 的 信 口 条 件 。 其 中 氧 的 入 
吕 温 度 低 于 其 拟 沸腾 温度 , 氧 处 于 液 相 状态 且 在 后 
续 组 合 过 程 中 会 发 生路 临界 转变 。 


T 表 1 燃料 与 氧化 剂 入 口 条 件 
- Tab.1 Inlet condition of fuel and oxygen 

参数 LOx CH, 
VitHR/ (g + s7') 143.1 44.4 
温度 /K 85 288 
密度 /(kg * m~?) 1178.3 42.1 
速度 /(m » s7!) 3.6 63.7 

„400 mm Adiabatic wall 
‘C i 50 mm 


Constant pressure outlet 


Al2 Mascotte 燃烧 器 结构 示意 图 


Fig.2 Schematic diagram of Mascotte combustor 


ae 
AEE, Se ESPN HERES HOV RESEDA IT Ft RRRA AO FHT g2 


图 3 计算 网 格 中 心 横 截 面 


Fig.3 Central plane of computation grid 


2 计算 结果 及 分 析 


2.1 流体 热 物 性 计算 结果 


为 了 验证 1.2 节 混 合 物 热 物性 计算 方法 的 准确 
性 ,采用 美国 国家 标准 与 技术 人 研究 院 (NIST) 的 数据 
作为 对 比 基 准 。 图 4 给 出 了 4 种 环境 压力 下 ,75% 
甲烷 与 25% 氧 混合 物 的 密度 与 定 压 比 热 随 温度 的 
变化 ,总 体 而 言 模型 的 准确 性 良好 。 
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(b) 定 压 比 执 
图 4 甲烷 / 氧 混 合 物 的 热 物 性 计算 结果 


Fig.4 Thermophysical properties calculation 


results of the mixture of methane and oxygen 
由 计算 结果 可 看 出 超 临 界 压力 下 流体 热 物 理性 
质 的 特点 。 其 中 ,流体 密度 随 温度 的 变化 呈现 连续 
分 布 ,不 存在 低压 下 由 相 变 引起 的 密度 间断 。 定 压 
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比 热 在 混合 物流 体 临界 温度 附近 会 出 现 峰 值 ,意味 
着 此 时 流体 大 量 吸 热 后 温度 的 变化 依然 很 小 。 定 压 
比 热 的 峰值 与 密度 的 大 梯度 区 域 相对 应 ,共同 描述 
了 流体 在 超 临 界 压力 下 由 类 液体 过 渡 到 类 气体 的 过 
程 。 此 过 程 与 亚 临 界 条 件 下 流体 的 沸腾 过 程 类 似 ， 
因此 也 被 称 为 拟 沸腾 ”。 值 得 一 提 的 是 ,密度 的 梯 
度 以 及 定 压 比 热 的 峰值 在 压力 略 高 于 混合 物 的 临界 
压力 时 达到 最 大 ,而 随 着 压力 的 升 高 会 同时 减 小 , 表 
明 “ 拟 沸腾 "现象 的 逐渐 减弱 。 


2.2 FPV 方法 火焰 面 数 据 库 的 建立 


将 真实 气体 模型 写 入 Flame Master 程序 ,用 于 
求解 火焰 面 方程 。 图 5 为 1.5 节 计 算 条 件 下 层 流 火 
焰 面 数据 库 中 最 大 温度 随 当量 标量 耗 散 率 变化 的 S 
型 曲线 。 对 于 同一 个 标量 耗 散 率 存 在 2 个 火焰 最 大 
E e e 


fat A AGE IP SC, PA EKHE) ~ (18) 
的 解 。 稳 定 燃烧 分 支 构成 了 SLE 方法 所 使 用 的 火焰 


据 库 ,FPYV 方法 则 同时 包含 稳定 燃烧 与 不 稳定 
婚 砷 分 支 ,因而 前 者 的 火焰 面 数据 库 是 后 者 的 子 集 。 
亿 侈 不 稳定 燃烧 状态 并 且 使 用 进度 变量 描述 化 学 反 
IEP BoE: FPV 方法 能 够 捕捉 汕 流 火焰 局 部 熄火 
= rears 


不 稳定 分 支 


图 5 液 氧 /甲烷 对 撞 火 焰 最 高 温度 随 标 量 耗 散 率 变化 曲线 


Fig.5 Variation of maximum temperature with scalar 


dissipation rate in LO,/CH, counter flame 

图 6 展示 了 与 $ 型 曲线 相应 的 各 火焰 面 解 的 温 
度 在 混合 分 数 空间 的 分 布 。 灰 色 虚 线 标示 了 当量 比 
混合 分 数 的 位 置 。 图 中 蓝 色 虚线 对 应 $ 型 曲线 的 转 
折 点 ,其 上 方 为 稳定 燃烧 分 支 , 下 方 为 不 稳定 燃烧 分 
支 。 最 下 方 曲线 表示 火焰 濒临 熄灭 的 状态 ,由 于 计 
算 条 件 为 高 压 纯 氧 燃 烧 , 炸 火 温度 可 延伸 至 约 600 K 
的 较 低 值 。 

为 了 反映 真实 气体 效应 的 影响 ,图 7 给 出 了 分 
别 使 用 真实 气体 模型 与 理想 气体 模型 计算 的 对 撞 火 
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焰 温 度 与 密度 分 布 , 横 坐标 采用 对 数 坐标 以 突出 液 氧 
一 侧 的 物性 差异 。 图 中 阴影 区 域 表 示 压 缩 因 子 Z <1 
的 大 致 范围 。 该 区 域内 混合 分 数 小 于 0. 01 ,流体 以 
低温 液 氧 为 主 且 真实 流体 效应 显著 。 对 比 可 见 使 用 
理想 气体 异型 的 计算 结果 出 现 很 大 误差 ,在 氧化 剂 边 
界 一 侧 真 实 气体 算得 密度 约 为 pw nc =1 178 kg/m’, 
而 理想 气体 算得 密度 约 为 pw ic =254 kg/m ,这样 的 
误差 会 造成 流 场 计算 中 的 严重 错误 。 在 阴影 区 域外 
流体 温度 为 常温 或 高 温 ,压缩 因子 趋 近 于 1 ,使 用 理想 
气体 模型 计算 的 密度 误差 和 火焰 温度 误差 并 不 多 。 
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图 6 不 同 燃烧 状态 的 火焰 温度 分 布 


Fig.6 Flame temperature distribution of 


different combustion state 
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图 7 采用 真实 气体 与 理想 气体 模型 的 火焰 面 方程 解 


Fig.7 Solution of flamelet equations using 


real-gas and ideal-gas model 
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男 一 方面 ,理想 气体 模型 没有 反映 出 低温 液 氧 
的 “ 拟 沸腾 "现象, 即 流体 从 类 液体 跨 过 临界 点 而 转 
变 为 类 气体 的 过 程 。 当 液 氧 达到 临界 温度 附近 时 定 
压 比 热 出 现 峰值 ,意味 着 此 时 即使 吸收 大 量 热量 , 温 
度 也 不 会 出 现 明显 变化 。 因 此 真实 气体 火焰 面 计 算 
结果 的 温度 比 理想 气体 火焰 面 更 低 , 差异 集中 于 混 
合 物 临界 温度 附近 ,最 高 温度 差 约 200 K, 


2.3 流 场 计算 结果 与 分 析 


2.3.1 液 氧 甲 烷 火 焰 的 典型 特征 

图 8 将 OH 的 数值 模拟 结果 与 实验 测量 的 OH 
辐射 的 Abel 变换 时 均 图 像 中 进行 了 对 比 。 可 见 数 
值 模拟 与 实验 结果 相 比 火焰 长 度 偏 长 ,相对 误差 为 
15% 左 右 。 此 外 ,模拟 结果 的 火焰 扩张 角度 偏 小 并 
目 友 应 区 域 的 宽度 也 更 窗 。 数 值 模拟 误差 的 原因 主 
要 伍 括 两 方面 :其 一 ,雷诺 平均 方法 无 法 对 强 剪 切 流 
APP RJ A EA 
ACES ELA WAS FE HE, REP Da SPL HE LEE TF 
(ERB 5 BABE WOKE A, SI TOE AE 
流 混合 预测 的 难度 。 总 体 来 看 ,仿真 结果 能 够 大 至 
准确 地 描述 超 临界 条 件 下 液 氧 甲烷 火焰 特有 的 锥 形 
SOM AYA. 
© m 


No 0.01 002 0.03 0.04 0.05 0.06 0.07 
上 da + + + he ate 


OH 
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图 8 数值 模拟 的 OH 质量 分 数 与 实验 的 
OH 辐射 时 均 图 像 对 比 

Fig.8 Comparison between the OH mass fraction of numerical 

simulation and averaged OH emission of experiment 

图 9 给 出 了 流 场 中 温度 .混合 分 数 的 分 布 情况 。 
由 混合 分 数 云 图 可 看 出 ,中 心 的 低温 液 氧 从 喷嘴 喷 
出 之 后 ,在 与 同 轴 高 速 甲 烷 气 流 的 剪 切 作 用 下 进行 
Bike. WANIE 拟 沸腾 ”过程 并 转化 为 气态 ， 
中 心 液 氧 核 变 细 并 最 终 完 全 溶解 。 气 态 氧 与 甲烷 掺 
混 形 成 恰当 混合 比 的 可 燃 混 气 , 随 后 完成 燃烧 ,形成 
了 包 囊 在 液 氧 柱 周 围 的 狭长 锥 形 火焰 。 在 流动 向 下 
游 发 展 的 过 程 中 , 剪 切 层 随 着 滑 流 挫 混 的 作用 而 不 
断 增 厚 ,混合 分 数 达到 当量 值 f, =0.2 附近 的 区 域 


= 
AEE, JET URIS HERES HY REBUN ALE I Fee BIE’ ERAT! 831 


不 断 变 宽 ,因而 燃烧 反应 区 域 相应 变 宽 ,火焰 呈现 出 
一 定 的 扩张 角度 。 

云图 中 标记 出 的 流 线 展 示 了 燃烧 室内 的 流动 形 
态 , 其 中 最 明显 的 特征 是 燃烧 室 头 部 与 火焰 下 游 存 
在 一 前 一 后 2 个 回流 区 。 前 方 回 流 区 的 产生 主要 是 
由 于 燃烧 室 壁面 的 约束 ,回流 的 高 温 已 燃气 体 有 助 
于 加 热 液 氧 核 并 且 稳 定 火 焰 。 后 方 回流 区 的 形成 则 
与 液 氧 的 “ 拟 沸腾 ”过程 相关 。 从 流 线 图 可 见 , 在 轴 
向 位 置 约 6 倍 射 流 直 径 处 ,流动 垂直 于 火焰 向 外 。 
这 一 方面 是 因为 液 氧 在 进行 “ 拟 沸腾 ”过程 时 迅速 
影 胀 , 另 一 方面 是 由 于 燃烧 而 导致 的 气体 膨胀 。 正 
是 由 于 液 氧 剧烈 膨胀 而 产生 了 治 径 向 向 外 的 流动 ， 
进而 形成 后 方 的 回流 区 。 


(a) 温度 


(b) 混合 分 数 
图 9 流 场 标量 云图 


Fig.9 Contour of scalar in flow field 

2.3.2 燃烧 室 压力 对 火焰 的 影响 

由 2.1 节 的 热 物 性 分 析 可 知 , 当 压力 升 高 时 低 
温 流体 的 “ 拟 沸腾 ”现象 变 弱 , 热 物性 的 改变 会 使 后 
续 的 燃烧 过 程 产 生 相 应 变化 。 本 市 研 究 燃烧 室 压 力 
对 液 氧 甲烷 火焰 的 影响 , 室 压 选取 6 MPa、10 MPa 以 
及 15 MPa, 燃 料 与 氧化 剂 入口 温 度 与 1.5 节 相 同 , 质 
量 流量 比 保持 mo = Men, =0. 31 不 变 。 

图 10 给 出 了 3 种 室 压 下 的 温度 场 计算 结果 。 可 
见 室 压 6 MPa 时 火焰 温度 分 布 与 实验 条 件 的 5.6 MPa 
基本 相同 ,在 轴 向 约 0.06 m 位 置 ,由 于 液 氧 的 膨胀 
火焰 突然 开始 扩张 ,在 火焰 下 游 有 回流 区 的 形成 。 
室 压 上 升 至 10 MPa 和 15 MPa 时 火焰 长 度 持续 增 
加 ,火焰 下 游 回 流 区 持续 缩小 直至 消失 。 一 方面 , 室 
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压 的 升 高 使 得 气态 甲烷 密度 增加 ,在 保持 流量 不 变 
的 情况 下 甲烷 入 口 速度 降低 , 减 小 了 外 圈 甲 烷 与 中 
心 液 氧 射流 的 速度 差 ,使 得 两 股 流动 之 间 的 庙 流 混 
合作 用 减弱 ,火焰 根部 的 燃烧 强度 相应 减弱 。 由 于 
算 例 设 置 为 富 燃料 燃烧 ,上 游 未 消耗 的 氧 会 在 下 游 
继续 参与 燃烧 ,火焰 长 度 相 应 增加 。 男 一 方面 , 流 线 
图 显示 了 室 压 15 MPa 时 火焰 下 洲 回 流 区 的 消失 。 
从 热力 学 角度 看 , 室 压 的 增加 使 得 流体 “ 拟 沸 腾 " 效 
使 得 液 氧 膨胀 产生 的 径 向 流动 强度 减弱 ,无 法 形成 
回流 。 
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10 不 同室 压 下 的 火焰 温度 云图 


Fig. 10 Flame temperature contour of 


different chamber pressure 

图 11 与 图 12 分 别 给 出 中 心 轴线 上 的 温度 与 密 
度 分 布 。 可 看 出 当 燃 烧 室 压力 升 高 , 轴 向 温度 分 布 
的 峰值 更 高 且 更 靠 后 。 火 焰 温 度 的 上 升 一 方面 是 由 
于 高 压 时 气态 甲烷 密度 更 大 ,有 助 于 增加 化 学 反应 
速率 ; 另 一 方面 则 是 由 于 火焰 变 长 之 后 , 放 热 区 向 下 
游 移动 。 密 度 分 布 曲线 显示 , 室 压 升 高 后 液 氧 入 口 
密度 略微 增加 ,2 种 工 况 密度 的 下 降 趋 势 一 开始 基 
本 相同 。 由 于 室 压 6 MPa 时 液 氧 的 “ 拟 沸腾 ”现象 明 
显 ,密度 在 0.06 m 位 置 急剧 下 降 ; 而 室 压 15 MPa 时 
由 于 “ 拟 沸腾 ”过 程 的 减弱 ,密度 变化 更 平缓 。 
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Fig. 12 Density distribution curves on central axial line 


3 4 论 


针对 超 临 界 压力 下 的 流动 与 燃烧 过 程 ,将 真实 
流体 模型 植 入 火焰 面 方程 计算 程序 ,开发 了 基于 火 
焰 面 进度 变量 模型 的 超 临 界 燃烧 的 均 相 求解 器 ,并 
对 液 氧 甲 烧火 焰 进 行 数值 研究 ,得 到 以 下 结论。 

1) 本 研究 建立 的 真实 流体 模型 能 在 广泛 的 热力 
学 状态 范围 内 准确 描述 流体 的 热 物 性 。 高 压条 件 下 
低温 流体 展现 出 “ 拟 沸腾 ”现象 ,其 热 物 性 对 于 一 维 
层 流 对 撞 火 焰 的 计算 结果 影响 显著 ,理想 气体 模型 
在 火焰 温度 与 液 氧 密度 的 计算 上 出 现 较 大 误差 。 

2) 本 人 研究 所 开发 的 均 相 求解 器 可 以 准确 地 预测 
超 临界 环境 下 液 氧 甲 烧火 焰 的 反应 区 形态 。 液 氧 在 
外 侧 甲 烷 气 流 的 剪 切 作用 下 完成 跨 临界 转变 与 掺 
混 , 使 得 火焰 形成 包围 于 液 氧 核 的 狭长 锥 形 反应 区 。 
由 于 液 氧 在 升温 过 程 中 迅速 膨胀 , 火焰 出 现 向 外 的 
突然 扩张 并 在 其 下 游 形成 回流 区 。 

3 ) 燃 烧 室 压力 升 高 时 液 氧 的 “ 拟 沸腾 ”效应 减 
弱 ,其 升温 过 程 中 的 膨胀 幅度 有 所 减 小 ,使 得 火焰 长 
度 增 加 且 扩 张 角 减 小 ,火焰 下 游 回流 区 消失 。 
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